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ABSTRACT 

The  nonlinear  stability  of  a  compressible  flow  past  a  wedge  is  investigated  in  the  hyper¬ 
sonic  limit.  The  analysis  follows  the  ideas  of  a  weakly  nonlinear  approach  as  first  detailed 
by  Smith  (1979).  Interest  is  focussed  on  Tollmien-Schlichting  waves  governed  by  a  triple 
deck  structure  and  it  is  found  that  the  attached  shock  can  profoundly  affect  the  stability 
characteristics  of  the  flow.  In  particular,  it  is  shown  that  nonlinearity  tends  to  have  a  sta¬ 
bilising  influence.  The  nonlinear  evolution  of  the  Tollmien-Schlichting  mode  is  described  in 
a  number  of  asymptotic  limits  which  were  first  identified  by  Cowley  and  Hall  (1990)  in  their 
linearised  account  of  the  current  problem. 
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§1.  Introduction. 

This  investigation  is  an  extension  of  the  work  of  Cowley  &:  Hall  (1990)  who  were 
concerned  with  the  linear  stability  of  hypersonic  flow  around  a  wedge  of  small  angle.  Here 
we  aim  to  examine  the  effects  of  including  nonlinearity  in  the  solutions  for  the  viscous 
modes  of  instability  considered  in  Cowley  &  Hall. 

Recently,  there  has  been  a  resurgence  in  research  of  hypersonic  flows,  motivated  by 
the  interest  in  developing  high  speed  aircraft.  In  compressible  flow  work  it  is  necessary  to 
specify  a  fui)jtion2il  relationship  between  the  temperature  and  the  viscosity  of  ti.e  fluid. 
Frequently  Chapman’s  law  is  used  in  which  the  viscosity  is  taken  to  be  directly  proportionsil 
to  the  temperature.  Although  this  law  can  be  rather  unrealistic  for  high  speed  flows, 
stability  analysis  using  this  relationship  is  often  more  straightforward  than  for  any  other 
commonly  used  law.  Cowley  &  Hall  (1990),  hereafter  referred  to  as  CH,  used  Chapman’s  in 
their  work  and  we  shall  follow  their  lead  here.  A  much  more  accurate  relation,  Sutherland’s 
law,  is  used  by  Blackaby,  Cowley  Sz  Hall  (1990)  in  their  investigation  of  the  instability  of 
hypersonic  flow  past  a  flat  plate.  They  obtained  smaller  growth  rates  for  the  vorticity 
mode  under  consideration  compared  to  those  when  Chapman’s  law  was  assumed.  In  this 
instance  we  know  that  Chapman’s  law  is  a  bad  approximation  for  the  vorticity  mode  since  it 
is  concentrated  in  the  layer  where  there  is  substantial  variation  in  the  ba^ic  temperature. 
However,  since  we  are  going  to  be  concerned  with  the  stability  of  Tollmien-Schlichting 
waves  which  by  their  nature  are  dominated  by  viscous  effects  clo'ie  to  the  wall  it  is  be  to 
hoped  that  the  results  obtained  here  using  Chapman’s  law  will  not  be  vastly  different  to 
those  found  should  a  more  accurate  law  be  used. 

At  hypersonic  speeds  real  gas  effects  will  be  important  but  we  do  not  take  them  into 
account  in  this  study.  These  have  been  considered  by  Fu,  Hall  Sz  Blackaby  (1990)  in 
their  investigation  of  the  instability  of  Gbrtler  vortices  at  hypersonic  speeds,  in  which  they 
also  assumed  Sutherland’s  law.  The  fluid  was  taken  to  be  an  ideal  gas  and  the  effects 
of  gas  dissociation  on  the  stability  of  the  modes  considered  were  obtained.  The  Gdrtler 
instability  at  hypersonic  speeds  has  also  been  investigated  by  Hall  &  Fu  (1989)  and  Dando 
&z  Seddougui  (1991). 

As  already  mentioned,  in  the  present  study  attention  is  focussed  on  (viscous)  Tollmien- 
Schlichting  waves.  CH  investigated  the  effect  of  an  attached  shock  on  the  viscous  and 
inviscid  modes  in  the  compressible  boundary  layer  on  a  thin  wedge.  For  a  discussion 
of  the  inviscid  modes  the  reader  is  referred  to  Mack  (1987).  Smith  &  Brown  (1990) 
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have  investigated  inviscid  inodes  at  high  Mach  numbers  and  in  the  absence  of  a  shock. 
These  authors  concentrated  on  the  vorticity  mode  while  CH,  in  their  study  of  the  inviscid 
modes,  looked  at  the  acoustic  modes  which  differ  from  the  vorticity  modes  in  as  much  as 
they  are  not  concentrated  within  a  thin  layer.  Blackaby,  Cowley  &  Hall  (1990)  discuss 
the  relationship  which  exists  between  these  various  modes.  For  large  Mach  numbers  the 
inviscid  modes  have  larger  growth  rates  them  the  viscous  ones  but  nevertheless  the  effect  of 
the  shock  on  the  vis.-ous  modes  is  of  some  importance.  In  fact  CH  show  that  the  presence 
of  the  shock  greatly  affects  the  viscous  modes. 

The  wavelengths  and  frequency  of  the  viscous  instability  in  a  compressible  boundary 
layer  at  hypersonic  speeds  were  given  by  Smith  (1989).  These  form  the  basis  of  the 
scalings  adopted  in  the  present  analysis.  CH  described  the  interactive  triple-deck  structure 
governing  disturbances  to  a  hypersonic  flow  over  a  wedge.  The  basic  flow  is  given  by  an 
exact  steady  solution  of  the  governing  equations  for  a  constant  density  gas  and  we  follow 
CH  and  assume  for  simplicity  that  the  wedge  is  insulated.  As  mentioned  by  CH  the 
zmalysis  can  be  extended  to  consider  an  isothermal  flow  where  heat  transfer  at  the  surface 
is  allow'"d.  CH  were  concerned  with  the  stability  of  the  basic  flow  to  viscous  modes  and  as 
a  first  analysis  they  sought  conditions  under  which  non-parallel  effects  could  be  neglected. 
They  showed  that  in  order  to  do  this  it  is  necessary  to  make  the  Newtonian  approximation 
and  this  yields  ?  “r  bound  on  the  Mach  number.  A  discussion  of  when  non-parallel 
effects  become  important  for  hypersonic  flow  is  given  by  Smith  (1989)  and  this  gives  the 
same  lower  bound  for  the  Mach  number.  CH  also  selected  the  angle  between  the  surface 
of  the  wedge  and  the  shock  so  that  *he  shock  occurs  in  the  upper  deck  of  the  triple-deck 
structure.  In  this  way  the  influence  of  the  shock  wave  on  the  Tollmien-Schlichting  waves 
can  be  investigated.  As  in  classical  triple-deck  analy.sir  as  the  perturbation  to  the  basic 
flow  is  increased  in  amplitude  from  an  infinitesimal  size  to  a  point  where  nonlinearity  must 
be  accounted  for,  it  is  the  equations  governing  the  disturbance  in  the  lower  deck  which 
are  the  first  to  become  nonlinear.  We  choose  to  retain  the  nonlinearities  in  the  lower  deck 
equations  since  the  aim  of  the  present  study  is  to  investigate  these  effects.  The  above 
requirements  give  bounds  on  the  position  of  the  shock  or  alternatively  bounds  on  the  size 
of  the  local  Mach  number  with  respect  to  the  size  of  the  Reynolds  number,  which  will  be 
taken  as  large  for  the  following  analysis. 
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CH  determined  the  linearised  conditions  which  hold  at  the  shock.  Since  the  nonlinear 
effects  are  unimportant  in  the  upper  deck  (where  the  shock  occurs)  to  the  order  of  approx¬ 
imation  considered  both  by  CH  and  the  present  study  these  linearised  shock  conditions 
are  appropriate.  It  is  found  that  the  effect  of  the  shock  is  weak. 

The  linear  results  of  CH  for  the  Tollmien-Schlichting  waves  will  be  presented  in  §4, 
however,  we  give  a  brief  summary  of  their  results  here.  CH  determined  the  nonlinear  system 
g  verning  the  instability  of  Tollmien-Schlichting  waves  over  a  thin  wedge  in  the  presence  of 
a  shock.  They  linearised  this  system  in  order  to  determine  the  dispersion  relation  relevant 
to  small  disturbances  and  neutral  solutions  were  presented  for  various  positions  of  the 
shock.  Since  the  disturbances  are  three-dimensional  there  exists  an  infinity  of  solutions. 
CH  presented  asymptotic  solutions  for  the  neutral  wavenumbers  which  agree  well  with 
their  numerical  results.  In  particular,  they  showed  that  as  the  position  of  the  shock,  y,, 
tends  to  infinity  the  shockless  dispersion  relation  derived  by  Smith  (1989)  is  recovered. 
CH  also  gave  a  discussion  on  the  temporal  growth  of  the  linearised  modes. 

CH  proceeded  to  investigate  the  effect  of  the  shock  on  neutral,  two-dimensional  invis- 
cid  modes  of  instability  having  wavelengths  scaled  on  the  boundary  layer  thickness.  These 
modes  are  concentrated  in  the  temperature  adjustment  layer  centred  at  the  generalised 
inflection  point  and  solutions  were  given  for  the  acoustic  modes  alluded  to  earlier.  (A 
discussion  of  the  vorticity  modes  in  the  absence  of  a  shock  is  given  by  Smith  &  Brown 
(1990).)  However,  CH  show  that  the  shock  has  a  negligible  effect  on  the  acoustic  modes 
unless  the  distance  of  the  shock  from  the  surface  is  of  the  order  of  the  boundary  layer 
thickness. 

Before  we  give  an  outline  of  the  paper  we  mention  some  additional  recent  investigations 
of  instability  of  hypersonic  flow.  Balsa  &  Goldstein  (1990)  investigi  ted  the  high  Mach 
number  instability  of  supersonic  mixing  layers  using  asymptotic  theory  They  use  the 
WKB  method  and  solve  the  compressible  Rayleigh  equation  for  a  aon-isothermal  free 
shear  flow,  so  in  a  sense  their  solutions  are  simileir  to  CH.  This  work  is  also  related  to  that 


on  the  vorticity  mode  by  Smith  &  Brown  (1990).  An  investigation  of  nonlinear  effects  — . . 

on  the  acoustic  modes  at  hypersonic  speeds  described  by  CH  is  given  by  Goldstein  & 

Wundrow  (1990).  Many  aspects  of  theoretical  research  into  hypersonic  flow  are  discussed  ^  i 

in  the  recent  paper  by  Brown  et  al.  (1991).  ^ — - — 

The  plan  of  the  paper  is  as  follows.  In  §2  we  describe  the  basic  flow  situation  for 
hypersonic  flow  around  a  thin  wedge,  and  give  the  equations  governing  the  flow.  The ' 


scalings  are  chosen  so  that  we  axe  able  to  neglect  non-parallelism,  to  include  the  effects  of 
nonlineeirity  in  the  lower  deck  of  the  triple-deck  structure  and  to  ensure  that  the  shock  is 
located  in  the  upper  deck.  The  triple-deck  structure  for  hypersonic  flow  is  described  in  §3 
and  the  scalings  chosen  to  eliminate  some  of  the  variables  for  convenience.  In  addition, 
we  state  the  linear  conditions  which  hold  at  the  shock  (and  which  are  described  in  detail 
in  CH).  In  ^4  we  consider  the  nonlinear  solutions  of  the  governing  system  of  equations. 
Rather  than  attempt  a  full  numerical  solution  of  the  nonlinear  equations  we  execute  a 
weakly  nonlinear  analysis  of  the  system  in  a  manner  virtually  identical  to  that  of  Smith 
(1979)  who  studied  the  nonlinear  stability  of  an  incompressible  boundary  layer  over  a 
flat  plate.  In  this  way  we  derive  an  amplitude  equation  of  classical  Stuart- Watson  type  to 
describe  the  evolution  of  the  disturbance.  Solutions  of  this  evolution  equation  are  discussed 
in  §5  where  we  also  draw  some  conclusions. 


§2,  The  basic  flow. 

We  consider  the  flow  of  a  compressible  fluid  over  a  wedge  of  semi-angle  0.  As  in  CH, 
we  suppose  that  the  supersonic  flow  has  velocity  magnitude  U  and  we  consider  the  situation 
when  the  wedge  is  symmetrically  aligned  with  the  oncoming  flow.  Then  symmetrical  shocks 
develop  on  either  side  of  the  wedge.  We  denote  the  semi-angle  of  the  shocks  by  cr  and  so 
the  angle  between  the  shock  and  the  wedge  is  ^  =  a  —  This  notation  is  chosen  to  be 
consistent  with  CH  and  the  situation  is  illustrated  in  Figure  1. 

We  take  {x,y,z)  co-ordinates  with  x  denoting  the  distance  along  the  upper  surface 
of  the  wedge,  y  the  distance  normal  to  the  wedge  face  and  z  as  the  spanwise  co-ordinate. 
We  suppose  that  the  corresponding  velocities  are  {u,  v,  w)  and  denote  quantities  upstream 
of  the  shock  and  in  the  region  between  the  shock  and  wedge  by  the  subscripts  u  and  s 
re-spectively. 

We  assume  that  the  fluid  is  a  perfect  gas  so  we  have  the  equation  of  state 

p  =  (7-l)^pf,  (2.1a) 

where  p,  p  and  T  denote  the  pressure,  density  and  temperature  of  the  fluid  respectively, 
whilst  7  is  the  ratio  of  specific  heats  and  Cp  is  the  specific  heat  at  constant  pressure.  Then 
the  speed  of  sound  in  the  upstream  flow  is  defined  by 

=  ^  =  1)K,  (2.1b) 

Pu 
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where  =  CpT  is  the  enthalpy  of  the  fluid.  It  follows  that  the  upstream  Mach  number  is 
given  by 

.  (2  Ic) 


The  inviscid  solution  for  hypersonic  flow  over  a  wedge  is  well  known  (see  for  example 
Hayes  &  Probstein  (1966))  for,  if  we  deflne  e  =  Pu|ps^  then  we  have 


(7  +  1 )  (7  -  sin^  (t)  ’ 


=  1  +  7M„(1  —  e)  sin^  cr, 


-p  =  1  +  i(7  -  1)(1  -  sin^  a, 
tan  <f>  =  €  tan  a, 


(2.2a) 


(2.26) 


(2.2c) 

(2.2d) 


(“1|)m  =  (^ll)a  =  («x)„  = -^sin<r,  (ujl),  = -ef^sina,  (2.2e) 

where  uy  and  u±  are  the  velocity  components  parallel  and  perpendicular  to  the  shock. 
Thus  we  see  from  (2.2e)  that  the  magnitude  of  the  fluid  velocity  in  the  region  between  the 
shock  and  the  wedge  is  given  by 

Ug  =  U  cos  <7  (1  +  tan^  cr)  *  ,  (2.3) 

which  yields  the  slip  velocity  along  the  wedge.  Thence,  from  (2.1c),  (2.2c)  and  (2.3)  the 
Mach  number  in  the  shock  layer  is  given  by 


2  2M^  cos^  O’  (1  +  tan^  a) 

"  ""2  +  (7-l)(l-e2)M2sin2(7’ 


whence 


Ml  = 


(2.4a) 


(2.46) 


“  2 cos^  <T  (1  +  tan^  ~  (7  —  ~  e^)M^  sin^  a 

Since  we  are  interested  in  Tollmien-Schlichting  instabilities  of  this  basic  flow  we  need 
to  consider  the  governing  equations  which  comprise  the  continuity,  Navier-Stokes  and 
energy  equations  for  a  viscous,  compressible  fluid.  We  non-dimensionalise  these  equations 
using  the  values  of  the  flow  quantities  between  the  shock  and  the  wedge.  We  write 

(x,y,2)  =  L(x,y,z),  iu,v,w)  =  Usiu,v,w), 
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where  the  lengthscale  L  is  the  distance  from  the  tip  of  the  wedge  to  the  position  of  interest. 
We  non-dimensionedise  time  with  respect  to  L/Ua^  pressure  with  respect  to  paUg  and  the 
other  variables  with  respect  to  their  values  in  the  shock  layer.  If  we  define  a  Reynolds 
number  for  the  fiow  by 

PsUaL 


Re  = 


where  pa  is  a  typical  viscosity  then  the  flow  is  governed  by  the  equations 

dp 


at 


+  V.(/>u)  =  0, 


Du  „  1 

P-7^  =  -^P  + 


2V.(,<«)  +  V  (iii'  -  |a-)V.u) 


Eind 


Dt  ^  '  Re  ^ 

=  (T,  _ 


(2.5) 

(2.6a) 

(2.66) 

(2.6c) 


where  p,  p  and  T  are  the  non-dimensional  density,  pressure  and  temperature.  The  shear 
and  bulk  viscosities  p  and  p'  have  both  been  non-dimensionalised  with  respect  to  pa  and 
Pr  is  the  Prandtl  number.  Finally,  the  components  of  the  rate  of  strain  tensor  are 

_  1  /  dui  duj\ 

2  \dxj  dxi )  ’ 

and  the  dissipation  function  $  is  given  by 

^  =  2pe:e+(^p'  -  (V.u)^  . 

In  the  following  analysis  we  shall  cissume  that  the  Reynolds  number  is  large.  From 
(2.1a)  the  non-dimensional  equation  of  state  is 


7-^aP  = 


(2.7a) 


and  the  non-dimensional  enthalpy  of  the  fiow  is  given  by 

h  =  T.  (2.76) 

As  in  CH  we  shall  restrict  attention  to  the  ceise  when  there  is  no  heat  transfer  on  the  wedge 
surface  so  that  dT/dy  =  0  when  y  =  0  and  the  Prandtl  is  taken  to  be  unity.  Then  the 
temperature  of  the  fluid  at  the  wall  is  given  by 

T„  =  l  +  i(T-l)M",  (2.8) 
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where  the  subscript  s  has  been  dropped  from  the  Mach  number. 

The  analysis  presented  in  this  paper  is  based  on  the  assumption  the  temperature 
dependence  of  the  viscosity  is  described  by  Chapman’s  law;  i.e. 

li  =  CT,  (2.9) 

where  C  is  a  constant.  The  argument  leading  to  the  choice  of  scalings  shown  below  was 
given  by  CH  for  the  more  realistic  Sutherland’s  law  as  well  as  for  Chapman’s  law.  However, 
since  the  results  presented  in  CH  are  restricted  to  the  Chapman  law  scenario  and  since  our 
attention  is  also  focussed  on  this  case  we  shall  just  give  the  scalings  appropriate  to  this 
situation  here. 

Then,  following  CH,  we  introduce  the  scalings 

cr  =  and  (7  —  (2.10a,  6) 

where  S  and  F  are  0(1)  constants.  If  (2.10a)  is  satisfied  for  a  small  then  the  shock  will 
occur  in  the  upper  deck  of  the  triple-deck  structure  and  the  effect  of  the  the  shock  on  the 
growth  rate  of  the  Tollmien-Schlichting  waves  can  easily  be  determined.  Consequently,  we 
are  assuming  that  <^,  the  angle  between  the  shock  and  the  wedge  face,  is  small. 

The  problem  is  greatly  simplified  if  the  contributions  from  non-parallel  effects  can  be 
neglected  in  a  rational  manner.  This  will  be  the  case  if  the  disturbance  develops  over  a 
smjill  streamwise  distance.  It  can  be  deduced  from  the  basic  flow  solution  (2.2a)  togef^her 
with  (2.4a)  that  e  <<  1  and  for  >  0  in  (2.46)  we  are  forced  to  make  the  ‘New'tonian’ 
assumption 

(7-l)«l,  (2.11) 

if  non-parallel  effects  are  to  be  ignored.  The  scalings  (2.10)  take  account  of  these  conditions 
and  these  considerations  furnish  a  lower  bound  for  a.  An  upper  bound  for  a  is  obtained 
by  including  the  effects  of  nonlinearity  in  the  lower  deck  which  is  the  aim  of  the  present 
study.  In  CH  the  linearised  problem  was  obtained  from  the  nonlinear  problem  they  initially 
formulated.  As  in  CH  we  intend  to  apply  linearised  shock  conditions  in  the  upper  deck  and 
this  consideration  whilst  simultaneously  demanding  nonlinearity  in  the  lower  deck  yields 
the  desired  upper  bound  on  a.  Specifically,  the  neglect  of  the  non-parallelism  and  the 
imposition  of  nonlinearity  imposes  the  constraints 

Re~'^  «  a  «  jRe"^,  (2.12a) 
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and  so,  from  (2.10a) 


Re'^  «  M  «  Re^.  (2.126) 

For  a  detailed  descriptive  of  the  argument  leading  to  (2.10)  and  (2.12)  the  reader  is  once 
more  referred  to  CH. 

If  the  upper  bou:  .J  on  cr  is  satisfied  then,  as  discussed  in  CH,  we  are  able  to  ignore 
viscous  effects  in  the  cniropy  and  shear  waves  which  are  inevitably  produced  when  the 

pressure/ acoustic  wave  meets  the  shock.  With  the  scales  described  above  and  assuming 

that  the  shock  has  a  viscous  internal  layer  the  entropy  and  shear  waves  will  have  y~  scales 
much  less  than  the  width  of  the  upper  deck. 

Prom  (2.106)  if  is  given  by  (2.8)  then 

Tu,  ~  irEV-2  »  1,  (2.13a) 

and  to  satisfy  >  0  in  (2.46)  we  require 

<  2.  (2.136) 

The  position  of  the  shock  is  given  by  (2.2d).  Thus,  in  cartesian  co-ordinates  at  the  shock 
we  have,  from  (2.2a)  and  (2.10a)  that 

y  —  xetancT  w  (2-14) 

We  re-emphasize  the  point  here  that  much  of  the  above  is  described  in  much  greater  deteiil 
by  CH. 


§3.  The  triple-deck  structure. 

It  was  shown  by  CH  that  the  hypersonic  flow  described  in  the  previous  section  is 
governed  by  a  triple-deck  structure.  In  deriving  the  relevant  scales  within  each  zone  it 
is  convenient  to  scale  out  a  number  of  the  variables.  The  triple  deck  structure  relating 
to  a  hypersonic  flow  situation  of  the  type  we  have  here  was  deduced  by  Smith  (1989)  in 
the  course  of  his  investigations  into  the  stability  of  supersonic  flows.  It  is  possible  to  use 
Smith's  results  to  infer  the  expansions  which  allow  /^u;,  and  M  to  be  eliminated  from  the 
governing  equations.  These  were  given  by  CH  but  will  be  repeated  here  for  convenience. 
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The  x—iZ—  and  t—  scalings  will  be  the  same  in  each  deck  and  are  given  by 

x  =  l  +  Re-»C^T^Mh-*X,  (3.1a) 


z  =  i?e-tctrjM-i7,  (3.16) 

and 

t  =  Re-^C^Tr^M^X-h,  (3.1c) 

where  C  =  fiwlTw  is  the  Chapman  constant  and  A  =  Ax“» ,  with  A  =  0.332,  is  the  Blasius 
boundary  layer  skin-friction  from  the  undisturbed  middle  deck  solution. 

§5.i  The  lower  deck 

The  lower  deck  is  the  region  in  which  the  nonlinearity  of  the  problem  manifests  itself 
and  the  equations  here  necessarily  contain  viscous  terms  since  it  is  this  layer  which  is 
required  to  reduce  the  disturbance  velocities  to  zero  at  the  wedge  surface.  Indeed  the 
relevant  equations  in  the  remaining  zones  of  the  triple-deck  are  linear  £is  we  shall  describe 
presently.  Within  the  lower  deck  we  have  scalings  which  take  the  forms 


(3.2a) 

Re-^C^T^M*\<U, 

(3.26) 

(3.2c) 

w  ^  Re-^C^T^M-^\hV, 

(3.2d) 

p  ~  +  Re-^C^M-^X^P, 

(3.2e) 

T  ~  Tu,, 

(3.2/) 

P  ~ 

We  substitute  these  expressions  into  (2.6a,  6)  which  give  at  leading  order 

(3.2y) 

Ux  +  Vy  +  Wz  =  0, 

(3.3a) 

Ur  +  UUx  +  VUy  +  WUz  =  Uyy, 

(3.36) 

II 

(3.3c) 
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Wr  +  UWx  +  VWy  +  WWz  =  -Pz  +  Wyy  •  (3.3c/) 

The  boundary  conditions  pertinent  to  system  (3.3)  require  no-slip  conditions  on  the  wedge 
face  so  that 

U  =  V  =  W  =  0  at  r  =  0,  (3.4a) 

and  the  solutions  must  match  with  those  in  the  middle  deck  as  Y  — >  oo  so  that 

U  —^Y  +  A(X,Z,r),  W—^0,  as  Y oo,  (3.46) 

where  .4  is  a  displacement  function  which  is  to  be  determined. 

^S.2  The  middle  deck 

The  middle  deck  covers  the  extent  of  the  undisturbed  boundary  layer.  Since,  from 
(2.13a),  Tfu  »  1  we  require  a  thin  transition  region  to  reduce  T  to  its  free-stream  value  of 
unity.  This  transition  region  has  thickness  O  ^(21ogru,)~^^  in  terms  of  the  Dorodnitsyn- 
Howarth  variable.  Thus,  the  middle  deck  is  itself  composed  of  three  regions:  (i)  &  bound¬ 
ary  layer  region  where  T  >>  1,  (it)  the  transition  region  and  (Hi)  the  region  above  the 
transition  region  where  T  1. 

The  thickness  of  region  (i)  is  O  (^Re~^C^Tw^  and  if  we  define 

y  =  Re-^C^T^y*,  (3.5) 

then  to  match  with  the  lower  deck  solutions  as  y*  — >  0  we  find  that  the  solutions  take 
the  forms 

«  ~  w5(y*)  +  Re-^C^TjMh-Uu^^.,  (3.6a) 

V  ~  -Re-^C*M-h^Axuoy  (3.66) 

w  ~  Re-^cKu-h^  (Po^oF^  (3.6c) 

(3.6c/) 

P  ~  Roiy*)  +  Re-^C^TSM'^X-^AR*oy».  (3.6.) 

Here  u^y*)  and  Ro(y*)  are  the  undisturbed  velocity  and  density  profiles  respectively  and 
Dx  =  —Pz-  For  Chapman’s  law  uj  takes  the  form  of  a  Blasius  function  for  a  compressible 
flow. 
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In  region  (ii)  we  write 


y  =  yc  -  i2e-^C^T„,(2logr«,)-^r*,  (3.7) 

where  yc  is  a  constant  which  determines  the  position  of  the  transition  layer.  Here  the 

solutions  are  very  similar  to  (3.6)  with  the  T^,  scales  slightly  altered  for  u  and  p.  In 
1  1 

particular,  for  u,  Tj  is  replaced  by  Tu,  ^  (21ogru,)^.  The  solutions  in  region  (in)  are 
similar  to  those  in  region  (i)  with  the  simplification  that  uq  =  Rq  =1  here. 

%S.3  The  upper  deck 

Finally,  in  the  upper  deck  the  basic  flow  quantities  take  their  free-stream  values. 
Thus,  as  in  CH,  the  pressure/acoustic  waves  have  scalings  given  by 


y  =  i?e-ictrjM-h-ty, 

(3.8a) 

u  ~  1  + 

(3.8b) 

(u,io)  ~  , 

(3.8c,  d) 

p  ~  +  Re~*C^M~^ X^p, 

(3.8e) 

p^l  +  Re-*C*M^X^p, 

T  ~  1  +  (7  -  l)i2e"5C^MU^f. 

(3.8y) 

If  we  substitute  (3.8)  into  (2.6)  and  (2.7)  and  rearrange  the  resulting  equations  we  find 
that  the  pressure  p  satisfies 

Pxx  -  Pyy  -  Pzz  =  0.  (3.9) 

The  solutions  in  the  upper  de  k  must  match  with  those  in  the  middle  deck  as  y  — y  0. 
From  the  y— momentum  equation  we  obtain  the  boundary  condition  that 

Py  =  Axx  at  y  =  0,  (3.10) 

while  p  —  P  when  y  =  0.  The  remaining  boundary  condition  is  that 

p  =  0  at  y  =  y,,  (3-11) 

where  y,  is  the  position  of  the  shock.  The  analysis  leading  to  this  condition  with  our 
basic  state  is  given  in  the  Appendix  of  CH  and  will  not  be  repeated  here.  The  reader  is 
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also  referred  to  Moore  (1954),  Ribner  (1954)  and  McKenzie  &  Westphal  (1968)  for  the 
conditions  which  hoH  at  the  shock.  This  condition  is  obtained  from  the  linear  inviscid 
problem  and  will  continue  to  hold  for  the  nonlinear  analysis  presented  in  this  paper  since 
the  nonlinearity  is  confined  to  the  lower  deck  to  the  orders  considered. 

For  Chapman’s  law  the  position  of  the  shock  is  given  by  (2.14).  Thus,  from  (3.8a), 
(2.13a)  and  (2.10a)  we  obtain 


els' 


(3.12) 


We  shall  now  follow  the  convention  of  referring  to  the  solutions  of  (3.9)  as  acoustic  waves. 


§4.  The  nonlinear  solution. 

In  general  the  solution  of  equations  (3.3),  (3.4)  and  (3.9-11)  is  a  fully  numerical  task 
but  it  is  possible  to  make  analytical  progress  by  implementing  a  weakly  nonlinear  analysis 
of  the  system.  The  method  used  here  is  identical  in  spirit  to  that  proposed  by  Smith 
(1979)  who  investigated  the  nonlinear  stability  of  an  incompressible  Blasius  boundary 
layer  to  Tollmien-Schlichting  type  disturbances.  We  differ  slightly  from  Smith  (1979) 
by  considering  a  three-dimensional  mode  as  in  CH  but  this  poses  no  formal  difficulties 
whatsoever. 

Our  objective  is  to  monitor  the  streamwise  development  of  a  finite  amplitude  TS  wave 
whose  leading  order  wavenumbers  in  the  streamwise  and  spanwise  directions  are  a  and  ^ 
respectively  and  whose  frequency  is  Consequently,  we  seek  modes  whose  fundamental 
component  is  proportional  to 

E  =  exp 

If  the  weakly  nonlinear  disturbance  is  allowed  to  develop  in  the  vicinity  of  the  linear  neutral 
point  and  if  the  relative  amplitude  of  the  disturbance  in  the  lower  deck  is  0(h),  h  -C  1, 
then  following  Smith  (1979)  it  is  found  that  the  amplitude  A  of  the  mode  will  evolve  on 
a  O(h^)  lengthscale.  However,  there  is  a  lower  bound  on  the  possible  size  of  h  which  is 
derived  by  consideration  of  the  non-parallelism  of  the  flow.  Using  an  argument  formally 
identical  to  that  of  Hall  Sz  Smith  (1984)  it  is  straightforward  to  demonstrate  that  the 
following  account  is  valid  for  all  h  in  the  range  0{Re~'^M'^T^  )<<h<<l.  If  the  lower 
bound  of  this  inequality  is  attained  then  the  lengthscale  over  which  the  amplitude  of  the 


(aX  +  fiZ  -  j\{q)dqj 


(4.1) 
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perturbation  modulates  is  identical  to  that  over  which  the  non-parallelism  of  the  basic  flow 
occurs.  For  full  details  of  this  aspect  the  reader  is  referred  to  Hall  &  Smith  (1984). 

Prom  (3.1a)  we  observe  that  the  linear  neutral  stability  point  has  been  non-dimensional- 
ised  so  as  to  be  at  a;  =  1  and  so  we  need  to  consider  perturbations  at  the  point 

x  =  l-f/i^a:2.  (4.2a) 

Since  the  skin  friction  A  is  a  function  of  x  it  too  will  be  slightly  perturbed  from  its  neutral 
value  Ai  according  to 

A  =  Ai(l  +  h^Aa).  (4.26) 

Additionally,  we  choose  to  fix  the  frequency  but  allow  the  spanwise  wavenumber  to  vary 
and  write 

0  =  (4.2c) 

These  perturbations  imply  that  for  <<  1  we  seek  solutions  of  the  lower  deck  equations 

(3.3-4)  and  the  upper  deck  problem  (3.9-11)  in  the  forms 

3 

{U,  V,  W,P,  A,P)  =  ((1  +  h'‘X,)Y,  0)  +  ^  (4.3) 

n=l 

Finally,  the  boundary  condition  (3.46)  which  describes  the  matching  with  the  middle  deck 
solutions  now  assumes  the  form 


a— »(i  +  /(“A2)(y+A),  IV— >0  as  r— >co. 


(4.4) 


To  account  for  the  slow  modulation  of  the  amplitude  on  streamwise  lengthscales  we 
introduce  the  co-ordinate  X  =  h^X  and  then  by  use  of  multiple  scales  we  replace  all  X 
derivatives  throughout  according  to  the  recipe 


d 

dX 


dX 


+  h^ 


d 

dX 


(4.5) 


Substitution  of  (4.3-5)  in  (3.3-4)  leads  to  a  hierarchy  of  problems  at  increasing  orders  in 
h.  We  address  these  problems  in  turn. 
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%4-i  The  first  order  'problem 


The  solution  for  the  0{h)  terms  in  (4.3)  has  the  linear  stability  form 

(UuVuWuPuAi,Pi)  =  (u,(X,Y)Mx,y),Wi(X,Y),Pi{X),Ai{X),p^{X,g))  E+c.c., 

(4.6) 

where  we  recall  that  E  is  defined  by  (4.1)  and  here,  as  in  the  remainder  of  the  paper,  c.c. 
denotes  complex  conjugate. 

The  problem  for  the  unknowns  in  (4.6)  is  precisely  the  linear  problem  considered  by 
CH.  It  is  convenient  to  define 


^  =  (fa)  3 


‘(-9 


and  then  we  obtain 


aVi  +  ft#,  =  ‘f/’*  r  Ai(q)dq, 

(^a)3A^'(£o)  Jio 


isQ, 

£o  14 

gs 

(4.7) 

Ai{q)dq, 

(4.8) 

where  Af(£)  denotes  the  Airy  function.  In  order  to  satisfy  the  boundary  conditions  as 
£  — ►  oo  then 

a5Af'(£o)-4i 


(4.9) 


with  K  =  Ai{q)dq.  From  the  continuity  equation  (3.3a)  we  deduce  that 


Vi  = 


(fa)3  Ai 


K 


Af'(£)  -  Af'(£o)  -  £  /  Ai{q)dq 


(4.10) 


Consideration  of  the  upper  deck  problem  (3.9-11)  yields  the  solution 
^2  sinh  -  a^ )5  (y,  -  y)j 


Pi  = 


and 


Pi 


\/W~^  cosh  -  a2)iy,j 

^2  sin  [(a^ -/?2)l(y,  _y)] 


■Ai ,  for  fii  > 


(4.11a) 


Ai,  for 


(4.116) 


■y/a2  -  fil  cos  j^(a2  - 

Note  that  the  solution  (4.116)  for  pi  does  not  decay  as  y^  — >  oo.  Since  pi  =  Pi  at  y  =  0 
we  obtain  from  (4.9)  and  (4.11)  the  dispersion  relation 

(fa)5/9^K  —  g^)^ 


tanh  j^(/?j  —  a^)^ya^ 


for  Pf  >  a. 


(4.12a) 
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for  /Sj  < 


(4.126) 


tan 

As  mentioned  by  CH  in  the  limit  y,  — ►  oo  for  /Jj  >  the  relation  (4.12a)  becomes  that 
obtained  by  Smith  (1989)  for  hypersonic  flow  with  no  shock. 

The  neutral  solution  of  (4.126)  for  a  and  /?i  real  is  determined  from  the  eigenrelation 

(I3l  -  a!^)i  =  020'^/?!  tanh  ^(/3?  -  Q!^)^y,)  ,  (4.13) 

and  ^0  =  — cii^  where  ci  «  2.3  and  C2  w  1.0.  The  corresponding  expression  for  /Jj  < 
is  similar  to  (4.13)  and  CH  obtained  neutral  solutions  to  these  two  equations  for  various 
values  of  yj.  The  dispersion  relation  admits  an  infinity  of  solutions  and  in  Figure  (2)  we 
present  some  of  these  solutions  for  the  three  cases  y*  =  1,  4  and  16.  For  each  y*  we 
observe  that  there  is  one  mode  for  which  a  — >  0  as  Pi  — y  oo  and  an  infinity  of  modes 
all  of  which  asymptote  to  a  =  Pi  as  Pi  — y  oo.  The  Mach  cone  is  defined  by  the  line 
ot  =  Pi  so  that  above  this  line  the  acoustic  waves  are  sinusoidal  and  supersonic  whereas 
below  this  line  the  modes  are  subsonic.  Hence  for  nearly  all  spanwise  wavenumbers  there 
are  an  infinity  of  supersonic  modes  and  precisely  one  subsonic  mode.:  We  shall  infer  precise 
details  concerning  the  asymptotic  limit  Pi  — y  oo  in  Section  5  where  we  shall  also  discuss 
solutions  of  the  amplitude  evolution  equation  for  i4i(X)  which  is  to  be  derived. 

§^.;S.  The  second  order  problem. 

By  examining  O(h^)  terms  in  (3.3)  and  (3.9)  we  conclude  that  the  solutions  of  these 
equations  at  this  order  take  the  forms 

U2  =  U2iE^  +  U22  +  (4.14) 


with  similar  expansions  for  1^2)  W2,  P2,  A2  and  p2.  Here  the  terms  U21,  U22,  1^21?  1^22}  W21 
and  W22  are  functions  of  X  and  F,  p2i  and  P22  are  functions  of  X  and  y  and  the  remaining 
terms  are  functions  of  X  alone. 

From  0{h‘^E^)  terms  in  (3.3)  we  find,  using  manipulations  similar  to  those  described 
by  Smith  (1979),  that 


ocU2i+PiW2i  =  '-^{Ai)\ia)-"^  j^^Fiq)  +  Anq)  f  Ai{t)dt  dq+B2(X)  Aii2h)dt, 

”  (4.15a) 
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which,  using  the  continuity  equation  (3.3a),  leads  to 


^21  =  -2i(ia)-s  (a?72i  +  dt 

j£o 


(4.156) 


In  (4.15)  we  have  made  the  definitions 


where 


R{2H)  =  -2-^  [2Ai(0^*"(0  +  Ai'{^o)Ai'{0]  • 


(4.16a) 


(4.166) 


Imposing  the  boundary  condition  (3.46)  as  ^  — >  oo  in  (4.15a)  gives 


AiaAai  =  ^  (ii)^m)~5  /°°  F  +  At'(0  t  Ai{t)dt  di  +  B2  H  Ai{2h)di.  (4.17) 

^  ^  J(o  [  Jio  J  J^o 

From  0{h‘^E^)  terms  in  the  momentum  equations  (3.36,  d)  when  evaluated  at  F  =  0  we 
obtain  the  relationship 


=  2ift*P21, 


(4.18) 


and  elimination  of  B2  between  (4.15a),  (4.17)  and  (4.18)  yields  a  first  connection  between 
P21  and  A21.  A  second  relation  between  these  quantities  is  obtained  by  studying  the  upper 
deck  problem  (3.9-11)  where,  on  considering  0(6^ J5^)  terms  we  find  that 


^_4(^2_q,2)P2i  =0, 


(4.19a) 


with  associated  boundary  conditions 


P21  =  0  at  y  =  Vs, 

=  -Aa^A2i  at  y  =  0. 
oy 


(4.196) 

(4.19c) 


Solution  of  (4.19),  when  combined  with  the  matching  requirement  with  the  main  deck 
solution  that  P21  =  -P21  at  y  =  0  yields  the  relation 


(4.20) 


for  /3j  >  a*  with  an  analogous  expression  when  /Jj  <  This  second  expression  for 
■P2i/-A21  enables  J?2(X)  to  be  found,  as  previously  described.  The  resulting  expressions 
are  complicated  and  so  we  shall  not  state  them  here. 

Turni.TLg  to  the  0(h^)  mean  flow  terms  C/22, 1^2,  W22,  ■P22,  A'i2  and  P22  in  the  governing 
equations  we  obtain  the  results  that 

^22=0, 


{ia)^ 


a^c/; 


22 


+  A 


(  9vl‘^ ,,  ^  dwi 


d(2  •  Q(2  J 

This  latter  equation  has  solution 


=  ((*«)^) 


(c)' 


Vi  a- 


(4.21a) 
+  c.c.  (4.216) 


aC.^2+Atr22  =  (iar8a^|— I*  /  dt  f  r(q)d,,  (4.22a) 

^  J^o  Joo 


where 


J^o 


no  =  ii 

Applying  the  boundary  condition  (3.46)  as  ^  — >  00  implies  that 


+  C.C. 


(4.226) 


A22  = 


iZi4|2  r<if  f  rm- 

1  ^  Joo 


(4.23) 


Finally,  the  upper  deck  problem  (3.9-11)  yields  P22  =  P22  —  0. 


§.^.5  The  third  order  problem. 

A  ^ 

At  this  order  we  determine  the  amplitude  equation  for  the  unknown  function  Ai{X). 
This  equation  arises  from  a  solvability  condition  on  the  0{h^E)  terms  in  the  systems  (3.3) 
and  (3.9-11),  solutions  of  which  are  sought  of  the  form 


C/3  =  EUzx  +  E^Uz2  +  +  C/34  +  E-^uif  +  E-^ui^^  +  E-^U^l’ , 

with  similar  expansions  for  ^3,^3,  A3, P3  and  p^.  The  desired  solvability  criterion  is 
obtained  from  consideration  of  the  coefficients  of  0{h^E)  terms  in  the  lower  deck  equations 
(3.3)  and  this  leads  to 

A 

^{olUzi  +  ^il^3i)  -  +  ^iW3i)  =  -(*«)~®-^  +  -  ^2W\,  (4.24a) 


1— 3  rrtc) 
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where 


-GU,X)  =(m-)-^(^  -  eo)^(a^i  +  PiWi)  +  {ia)Hi  -  6)A2(aJ7i  +  PiWi) 

oX 

+  aA2Fi  +i^i^2Pi  +i(Q;^i  ■{•  ^\W\)^^\oiU2\  +/?iW^2i) 

1  d 

+  i{aU22  +  ^iW22){oiUi  +  +  {ioi)^Vi-^{aU22  +  ^iW^22) 

+  (ia)iFi''>^(aa2i  +  *1^2,)+ ((ia)i)‘'V2,^(a£j<'>  + 


(4.24i) 


The  boundary  conditions  obtained  from  (3.4)  are 


U3ii^o,X)  =  W3iiio,X)  =  0, 


(4.25a) 


Un{oo,X)  =  A3i(X)  +  A2ii(X),  P^3i(oo,X)  =  0. 


(4.256) 


By  evaluating  the  streamwise  and  spanwise  momentum  equations  (3.36,  d)  at  ^  =  fo  we 
get 

+iia-SAftA  =  (“^^  •  (^-26) 

The  homogeneous  forms  of  equation  (4.24)  and  boundary  conditions  (4.25)  are  sat¬ 
isfied  by  aUi  -F  fiiWi  and  so  (4.24-5)  only  has  a  solution  for  aUzi  +  ^iWzi  if  a  certain 
compatibility  condition  holds  (as  in  Stuart  (1960)  and  Watson  (I960)).  This  condition  is 
derived  by  considering  the  adjoint  system  of  the  problem  as  in  Hall  &  Smith  (1982).  We 
multiply  (4.24a)  by  the  function 


where 


K{o = Ai{o  - 


(4.27a) 


(4.276) 


and  integrate  over  the  range  (6,oo).  This  gives 


Az'(^o) 

C'Uo) 


aX2A\  +  (ia)  » 


K'(i)a((„x)di+ m)  di. 


(4.28) 


In  the  upper  deck  the  problem  for  pai  takes  the  form  (on  equating  terms  of  0{h^E) 
in  (3.9)), 


Qg2  "  (^1  “2)P3l  =  -t-2^l/32Pl, 

(4.29a) 

with  boundary  conditions 

p3i  =0  at 

1;% 

II 

(4.296) 

=  —a^Asi  +  2ia:^^  and 

dy  ab¬ 

P31  =  P31  at  y  =  0. 

(4.29c) 

solution  of  (‘’.•.29)  yields  a  second  expression  for 

i^a~iAi(^o)/3iP3i  + 

in  terms  of  the  solutions  derived  when  considering  the  first  and  second  order  problems 
in  Sections  4.1  and  4.2.  Then,  after  considerable  manipulation,  substitution  of  (4.245) 
for  G(^,X)  in  (4.28)  enables  the  desired  evolution  equation  to  be  obtained.  Since  the 
expressions  occurring  in  this  equation  are  long  and  complicated  we  have  relegated  these 
details  to  the  Appendix.  In  the  following  section  we  analyse  some  solutions  of  the  amplitude 
equation  and  comment  upon  some  of  the  implications. 


§5.  Results  and  discussion. 


By  implementing  the  procedures  described  in  the  previous  section  we  obtain  a  classical 

A 

evolution  equation  for  the  amplitude  Ai(A’)  which  takes  the  general  form 


dAi 


dX 


^  =  (ai^2  +a2^2)Ai  +a3Ai|Ail  , 


(5.1) 


where  the  definitions  of  the  complex  valued  constants  oi ,  02  and  03  are  given  in  the  Ap¬ 
pendix.  As  is  usual  in  this  type  of  v/eakly  nonlinear  analysis  interest  is  primarily  centred 
upon  the  sign  of  the  real  part  of  03.  It  is  easily  deduced  from  (5.1)  that 


_d 

dX 


^  =  [2Re(a2)^2  +  o,iP2)] \Ai\^  +  2Re{a3)\Ai\* ^ 


(5.2) 


and  thence  if  Re{a3)  >  0  nonlinear  effects  are  destabilising  and  the  equilibrium  state  with 
amplitude 


-Re(az) 


[ilfc(ai^2  +  02  ■^2)] 


(5.3) 
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where  n  =  1, 2, . . .  and  where  we  recall  that  C2  «  1.0.  Then,  using  the  functions  defined  in 
the  Appendix  it  follows  that  at  leading  orders  in  the  evolution  equation  for  Ai  becomes 


dX  yiAi{(o)  \(n  +  \)'K  j 


where  the  constants  Lj  are  defined  in  the  Appendix.  On  evaluation  this  amplitude  equation 
becomes 


- (-0.378,0.273)  +  v— - - 

y)(n+i)i  yl(n  +  i)t 


yl 


(-0.001,  -0.932)  +  — ^Vr(0, 0.318) 


(n  +  i)^ 


+  ^(-0.279,7.03  X  10"2)ii|ii!2. 

Va 


(5.46) 

We  may  deduce  from  this  equation  that  for  <<  1  all  the  coefficients  in  this  amplitude 
equation  become  small.  Further,  on  considering  the  coefficient  of  the  linear  term  we  see 
that  variations  in  the  perturbed  spanwise  wavenumber  P2  have  little  effect  on  the  linear 
growth  rate  of  the  TS  wave.  Finally,  for  ^02  =  0  we  have  the  possibility  of  a  finite  amplitude 
stable  mode  of  size 


(n  +  i)5 


A  second  case  for  investigation  is  that  of  supersonic  modes  with  a  >>  1.  From  Figure 
(2)  we  observed  that  such  neutral  modes  have  and,  more  accurately,  we  find  that 


-  f  nTrN  1  f  n'K\  1,1/  n7r\  .. 

^  V  ya  /  2A  V  y,  /  SPl  C2y,  V  y,  j  ^ 


10 

'  _j. 


(^5a) 


where  n  =  1,2, _  Thus  for  /?i  >>  1  the  possible  a  values  for  neutrally  stable  linear  TS 

waves  are  separated  by  0(l/n).  In  this  case  the  amplitude  equation  asymptotes  to 


PlylAi{^o)  dAi 
(n7r)2  die 


Az/Sf  L2  +  ^2 


s  (n7r)2 


(5.56) 
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Performing  the  necessary  evaluations  yields 


dAi 

dX 


(-8.016,5.781) 


n^/?i 


-¥ 


yl 


-Aj  +  ild2 


Ai  +  (-6.058,  -0.377)^^ 

y» 


(5.5c) 


and  hence  we  can  obtain  supercritically  stable  modes  of  amplitude  «s  1.15/?f  *\/— A2.  We 
note  that  in  this  a,  >  1  limit  the  scaled  equilibrium  amplitude  is  small  and  is  independent 
of  the  perturbed  spanwise  wavenumber  ^2. 

A  final  regime  for  asymptotic  investigation  is  that  of  the  solitary  subsonic  mode  in 
Figures  (2),  i.e.  those  modes  for  which  a  — >  0  as  — *  00.  The  solution  of  eigenrelation 
(4.12a)  in  this  case  assumes  the  form 


Cl 


(5.6a) 


and  then  the  amplitude  equation  (5.1)  simplifies,  at  leading  order  in  ^  to 


C2L\  —  2Ai(^o)  dA\ 


dX 


4* 


/?2(c2lf3  +  ^^4) 


-4i  + 


(02X5  {-2X6) 

^2^1  (cji?  +  SLg) 


aMiW 


or. 


^  =  [(-0.376, 0.117)A2i3f^  +  (0.125, -0.373)/92/Sf^] 

+  (-0.321,  -8.83  X  10"2)^-5ij \Ax  p. 


(5.66) 


This  evolution  equation  admits  a  stable  finite  equilibrium  solution  in  which  ..4i  =  0{l3i). 

We  remark  that  the  three  asymptotic  limits  briefly  alluded  to  above  agree  well  with  the 
numerical  solutions  sketched  in  Figures  (3)-(5).  It  is  straightforward  to  make  a  comparison 
of  the  relative  likelihoods  of  the  equilibrium  modes  in  the  »  1  case.  We  have  seen  that 
for  >>  1  the  supersonic  modes  have  amplitudes  O(/0i  ®)  whereeis  the  corresponding 
amplitude  for  the  subsonic  mode  is  much  larger,  0(/?i).  Hence,  in  practice,  at  large 
spahwise  wavenumbers  /?i  the  subsonic  mode  will  likely  dominate  the  other  modes  and  is 
therefore  much  more  likely  to  be  observable. 

To  conclude,  we  have  presented  an  account  of  some  properties  of  classical  weakly 
nonlinear  stability  theory  as  applied  to  an,  admittedly  idealised,  hypersonic  flow  problem 
originally  considered  by  CH.  Our  main  conclusion  is  that  for  all  the  parameter  regimes 
investigated  near  linear  neutrally  stable  modes  appear  to  be  supercritically  stable;  i.e.  the 
nonlinear  effects  tend  to  stabilise  the  slightly  unstable  linearised  modes. 
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There  is  much  further  work  which  is  possible  in  relation  to  this  problem.  The  most 
obvious  is  that  of  a  full  nonlinear  solutions  of  the  lower  deck  equations  (3.3)  together 
with  boundary  conditions  (3.4);  work  which  when  combined  with  the  relevant  upper  deck 
problem  clearly  needs  to  be  numerical  in  character.  A  second  course  of  action  is  that  of 
considering  viscosity/  temperature  dependences  other  than  that  of  Chapman’s  law.  As 
mentioned  earlier  CH  and  Blackaby  et,  al.  (1990)  consider  situations  in  which  the  more 
physically  realistic  Sutherland’s  law  is  used  although  CH  did  not  calculate  any  results  for 
this  case.  The  effect  of  changing  the  viscosity  law  on  our  results  would  be  of  great  interest 
for,  of  course,  if  it  transpires  that  the  properties  of  our  solutions  are  quite  sensitive  to  the 
actual  laws  invoked  then  this  raises  the  question  as  to  which  to  use  in  future  studies. 

One  complete  subject  area  which  has  been  deliberately  excluded  in  both  CH  and  the 
present  study  is  that  as  to  the  nature  of  any  nonlinearities  which  might  occur  in  the  vicinity 
of  the  shock  layer  position  y,.  The  scalings  chosen  here  were  such  that  (2.12a)  holds; 
technically  once  <t  reaches  a  value  of  0{Re~^)  then  the  shock  layer  problem  becomes 
nonlinear  itself  for  we  can  no  longer  ignore  viscous  effects  in  entropy  and  shear  waves 
which  are  produced  when  the  acoustic  wave  meets  the  shock.  Finally,  we  mention  the 
observation  of  CH  who  concluded  that  in  addition  to  the  neutral  modes  sketched  in  Figure 
(2)  there  are  other  modes  with  zero  growth  rate  but  with  infinite  frequency  and  therefore 
these  modes  cannot  be  described  within  the  asymptotic  framework  used  here.  Investigation 
of  both  the  linearised  and  nonlinear  problems  which  axise  within  the  necessarily  modified 
structure  would  be  extremely  worthwhile. 

There  are,  therefore,  many  directions  in  which  the  results  presented  here  may  be 
extended.  Progress  in  any  of  these  would  be  valuable  indeed. 
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APPENDIX 


In  this  appendix  we  summarise  the  results  of  manipulating  (4.28)  and  using  the  first 
and  second  order  solutions  in  order  to  derive  the  desired  amplitude  equation.  We  recall 
that  aUi  +  /?i  Wi  is  given  by  (4.8)  and  by  considering  the  0(h)  terms  in  the  lower  deck 
equation  (Z.Zd)  we  obtain  Wi.  Similarly,  W21  can  be  found  by  examining  0{h‘^E^)  terms 
in  this  equation  and  then  1/21  derived  via  (4.15a).  Finally,  Vi  and  V21  are  given  by  (4.10) 
and  (4.156)  respectively  and  thus  G(^,X)  deduced  from  its  definition  (4.246).  We  may  now 
proceed  to  substitute  in  (4.28)  and  the  resulting  classical  Stuart-Watson  type  evolution 
equation  for  Ai(a  )  is  as  follows: 

If  then 


+  a* 


r  /  L5  +  $+1(6 


2/7  +  $+2/8 


Ai\Ai\\ 


(A.la) 


where 


^  q2  -oc^Vs), 

and  where  {2/j}j._j  are  complex  constants  defined  presently.  If  however  then 


+  a5 


/ 2/5  +  $-2/6  ^  I  I  7  ,2 

(iTTwnrJ  ' 


(A.U) 


where 


24 


The  constants  here  are  defined  by 


{A.2a) 


,  i(ri)  2i*A!Ko)v4i”(6) 

■^2  = - > 

K  K 

L,  =  L({„)  -  rK(()Am  (  f  /’ di  , 

^  *'$0  \*'^0  / 

L4  =  -iA2(6)> 

^  (r20)(T14)  +  (r21)(T7) 

“  /c|k|2 

,  (r20)(ri3)  -  (T2i)(ri2) 

•^6  —  I.  10  •> 


{A.2b) 

(A.2c) 

(A.2d) 

(A.2e) 

(A.2/) 


i7  =  -2lf;^(^i(0)]  ,  (A,2«) 

and 

i*00 

Ls=2U'^  /  Ai(3)ds.  (A.2/1) 

J2t(o 

The  complex  constants  Tl,r7,ri2,ri3,T14,r20  and  r21  are  defined  and  evaluated 
in  the  Appendix  of  Bassom  (1989).  All  the  above  expressions  were  found  numerically. 
Initial  conditions  at  a  suitably  large  value  of  say  ^ooi  were  found  using  asymptotic 
methods  and  the  defining  equations  integrated  using  a  fourth  order  Runge-Kutta  scheme. 
The  integrations  necessary  in  the  evaluation  of  (A2)  were  performed  using  the  trapezium 
rule  combined  with  Richardson  extrapolation.  Various  values  of  ^00  and  step  length  were 
used  and  the  results  quoted  below  are  believed  to  be  accurate  to  within  0.1%.  Further 
checks  were  made  on  these  computed  values  by  using  the  program  to  calculate  some  of  the 
coefficients  given  by  Hall  &  Smith  (1982)  and  defined  in  their  equation  (2  27).  Using  the 
above  techniques  we  found  that 


Li  w  (0.4766,0.9858),  L2  w  (0.0277,-1.202), 
Lz  «  (1.916, 1.450),  I4  «  (-0.9573,  -0.7247), 
Ls  «  (-0.8346,1.389),  Le  w  (2.092,  -0.9825), 
I7  «  (-4.633,-0.6513),  Tg  «  (0.3473,1.536). 
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Figure  (1).  The 


geometry  of  the  wedge  and  shock  for  a  high  Mach  number  flow. 
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Figure  (2).  (a)  The  first  three  neutral  curves  a  =  a(y3i)  given  by  (4.12)  for  y,  =  1.  (b) 
The  first  seven  neutral  curves  for  ya=^‘  (c)  The  first  ten  neutral  curves  for  =  16. 
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Figure  (3).  The  real  part  of  the  coefficient  as  defined  in  (5.1),  (a)  The  values 

tahen  on  the  three  curves  in  Figure  (2a)  where  y*  =  1.  (b)  The  values  taken  on  the  seven 
curves  in  Figure  (26)  where  y,  =  4.  (c)  The  values  taken  on  the  ten  curves  in  Figure  (2c) 
where  y,  =  16. 


Figure  (3).  The  real  part  of  the  coefficient  ai(^i)  as  defined  in  (5.1).  (a)  The  values 
taken  on  the  three  curves  in  Figure  (2a)  where  y,  =  1.  (h)  The  values  taken  on  the  seven 
curves  in  Figure  (26)  where  ya  =  4.  (c)  The  values  taken  on  the  ten  curves  in  Figure  (2c) 
where  y,  =  16. 
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Figure  (4),  The  real  part  of  the  coefficient  02(^1)  as  defined  in  (5.1).  (a)  The  values 
taken  on  the  three  curves  in  Figure  (2a)  where  =  1.  (b)  The  values  taken  on  the  seven 
curves  in  Figure  (26)  where  y,  =  4,  (c)  The  values  taken  on  the  ten  curves  in  Figure  (2c) 
where  y«  =  16. 


Figure  (4).  The  real  part  of  the  coefficient  a2{0i)  as  defined  in  (5.1).  (a)  The  values 
taken  on  the  three  curves  in  Figure  (2a)  where  ys  =  1.  (b)  The  values  taken  on  the  seven 
curves  in  Figure  (26)  where  =  i.  (c)  The  values  taken  on  the  ten  curves  in  Figure  (2c) 
where  =  16. 


Figure  (4).  The  real  part  of  the  coefficient  a2(/?i)  as  defined  in  (5.1).  (a)  The  values 
taken  on  the  three  curves  in  Figure  (2a)  where  y,  =  1.  (b)  The  values  taken  on  the  seven 
curves  in  Figure  (26)  where  y,  =  4.  (c)  The  values  taken  on  the  ten  curves  in  Figure  (2c) 
where  y,  =  16. 
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Figure  (5).  The  real  part  of  the  nonlinear  coefficient  az{fix)  as  defined  in  (5.1).  (a)  The 
values  taken  on  the  three  curves  in  Figure  (2a)  where  =  1.  (b)  The  values  taken  on  the 
seven  curves  in  Figure  (26)  where  y,  =  4.  (c)  The  values  taken  on  the  ten  curves  in  Figure 
(2c)  where  y,  =  16, 
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is  shown  that  nonlinearity  tends  to  have  a  stabilising  influence.  The  nonlinear 
evolution  of  the  Tollmien-Schlichting  mode  is  described  in  a  number  of  asymptotic 
limits  which  were  first  identified  by  Cowley  and  Hall  (1990)  in  their  linearised 
account  of  the  current  problem. 
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